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COMPARISON OF STABILITY AND CONTROL PARAMETERS 
FOR A LIGHT, SINGLE-ENGINE, HIGH-WINGED AIRCRAFT 
USING DIFFERENT FLIGHT TEST AND PARAMETER ESTIMATION TECHNIQUES 


William T. Suit and Robert L, Cannaday 
SUMMARY 


t 

Longitudinal and lateral stability and control parameters were estimated 
from flight data for a high-wing, general aviation, airplane using flight data 
obtained at various flight conditions within the normal range of the aircraft. 
These parameters were estimated using an output error technique (maximum 
likelihood) and an equation error technique (linear regression). Longitudinal 
static parameters were also estimated from climbing, descending, and quasi- 
steady-state flight data. For the lateral excitations, four input forms were 
used involving some combination of rudder and ailerons. The resulting longi- 
tudinal and lateral parameter estimates were used to compute the periods and 
time-to-damp to one-half amplitude of the various aircraft modes of motion to 
determine the sensitivity of these motions to variations in the parameter 
estimates . 


INTRODUCTION 


The use of simulators for research investigations into aircraft dynamics 
is becoming an increasingly important tool of the research engineer. With 
this increased use of simulations comes the demand for greater simulator 
fidelity, which requires improved mathematical models. The aerodynamics of an 
aircraft can be described mathematically using its stability and control 
parameters in a set of equations of motion. 

Several techniques have been used in the past to estimate stability and 
control parameters from flight data including analog matching (refs, 1 and 2), 
the time vector method (ref, 1), and regression analysis (refs, 3 and 4), In 
recent years the estimation of parameters from flight data by use of a maximum- 
likelihood algorithm (ref. 5) has become fairly routine where adequate computer 
facilities are available. In general, the analyses reported in these 
references were based on the small perturbation equations of motion. 

To date, stability and control parameters for several low-winged light 
airplanes have been determined from both flight test data and wind-tunnel tests 
(refs, 6, 7, 8, and 9), For high-winged configurations wind-tunnel test data 
are available (ref, 10) but little has been published on the estimated values 
of stability and control parameters determined from flight data. Wind tunnel 
tests usually do not include estimated values for rotary derivatives, which can 

be estimated from flight test 'data using current estimation techniques, 
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The present paper is similar in some respects to the work that was done 
for the low-winged, general aviation aircraft as reported in references 7 and 8, 
However, in the present study, flight test data were obtained at three different 
airspeeds corresponding approximately to landing, approach, and cruise 
conditions. During the flight test, the aircraft was perturbed from trim 
conditions using either elevator or rudder and ailerons, and the stability and 
control parameters were estimated from the resulting data using both an output 
error method (maximum likelihood of ref. 5), and an equation error method 
(linear regression of refs. 3 and 11). Static longitudinal aerodynamic 
parameters were also calculated from steady climb, descent, and quasi-steady 
flight data for comparison with those obtained from perturbation flight data. 

This report describes the flight test procedure, presents the results 
obtained from the perturbation flight tests using the two estimation procedures, 
and compares values for the longitudinal static parameters with values 
calculated from independent flight tests at two different center-of-gravity 
locations. 


SYMBOLS 


The aerodynamic parameters are referenced to a system of body axes with 
the origin at the airplane center of gravity, which is located at 28.8 percent 
c, and with orientation .of body axes as shown in figure 1, which also shows 
the direction of positive forces, moments, displacements, angles, and linear 
and angular velocities. 
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Mx. My, M^ 

m 


acceleration measured along X, Y, and Z body axes, 
respectively, g units 

wing span, m 

wing mean geometric chord, m 

force along X, Y, and Z body axes, respectively, N 

acceleration due to gravity, m/sec^ 

moment of inertia about X, Y, and Z body axes, respectively, 
kg-m^ 

product of inertia, kg-m^ 

distance from airplane center of gravity to center of 
pressure of horizontal tail, m 

rolling, pitching, and yawing moments, respectively, N-m 

mass, kg 
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p 


U, V, w 


u' , v' , w' 


X, Y, Z 


X, y, z 


roll rate, rad/sec 
pitch rate, rad/sec 

2 

dynamic pressure, N/m 
yaw rate, rad/sec 

estimate of error covariance matrix 
wing area, m^ 

thrust, N -- f 

velocity along X, Y, and Z body axes, respectively, m/sec 

velocity component along X, Y, and Z body axes, 

respectively, at angle-of-attack sensor on wing-tip boom, 
m/sec 

control vector 

airplane total velocity, m/sec 

matrix, of measured states and input variables 

body coordinate axes through airplane center of gravity 

state vector 

output vector 

X-, y-, and z-coordinates, respectively, of the sensors on 
wing-tip boom relative to airplane center of gravity, m 

measurement vector 

angle of attack, rad 

angle of sideslip, rad 

left aileron deflection minus right aileron deflection, rad 

stabilator deflection, positive trailing edge down, rad 

rudder deflection, positive trailing edge left, rad 

angle between thrust axis and airplane X body axis, positive 
for thrust up, rad 

pitch angle, rad 
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tn 
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parameter vector 

o 

air density, kg/m 
measurement noise vector 
roll angle, rad 

perturbation in parameter vector 
lift coefficient, gm/qS 
rolling -moment coefficient, Mj^/qSb 
pitching-moment coefficient, My/qSc 
yawing-moment coefficient, M^/qSc 

thrust coefficient, T/qS (used in some publications as T^) 
axial-force coefficient, 
side-force coefficient, F^/qS 
normal -Jorce coefficient, F^/qS 
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Subscripts : 

c 

k 

m 

o 

t 


computed 

index 

measured 

coefficient at trimmed conditions 
trimmed conditions 


Superscripts : 


-1 

T 

M 

o 

A 


inverse matrix 
transpose matrix 
measured quantity 
nominal evaluation 
estimated value 

A dot over a symbol signifies a derivative with respect to time. 


DESCRIPTION OF AIRPLANE AND DATA SYSTEM 


The subject airplane was a four-place, externally braced high-wing, fixed 
tricycle landing gear, single-engine airplane, as shown in figure 2. Its 
pertinent geometric details and mass characteristics are given in table I. The 
mass characteristics shown were obtained from manufacturer's data on the subject 
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aircraft. The airplane instrumentation used to. record control-surface movements 
and airplane responses to these movements was basically like that of the subject 
aircraft of reference 7, 

The instrumentation system measured and recorded on tape the data used in 
this study. The variables recorded and the range of each sensing instrument 
is given in table II, The accuracy of these measurements is considered to be 
2 or 3 percent of full scale on each instrument. Unlike the data recorded in 
reference 7 in which FM and PAM were merged, all the data were recorded on FM 
(continuous) channels. The advantage of this was that all data channels could 
be filtered using an analog filter without introducing time delays in some 
channels relative to others. These data were digitized, then sampled at 
20 points per second and converted to engineering units to obtain the data used 
in this study. The pitot-static head for measuring velocity (dynamic pressure), 
and the angle-of-attack and angle-of-sideslip vanes were mounted on a boom 
located near the left wing tip. The boom extended 3/4 chord ahead of the wing 
leading edge. 


FLIGHT TESTS 


Three types of flight tests were flown to obtain the data used in this 
report, 

(1) Perturbation tests ,- These tests consisted of trimming the airplane 
with power for level flight, idle, or full power and perturbing the trimmed 
condition with either elevator or rudder and aileron doublets. A typical 
time history of the input forms used to excite the longitudinal motions is 
presented in figure 3. The aircraft was trimmed at three indicated airspeeds; 
31.6 m/sec (61 knots), 40,7 m/sec (78 knots), and 54,2 m/sec (104 knots), which 
correspond roughly to landing, approach, and cruise, although no flaps were 
used in these tests. 

The lateral perturbation tests consisted of trimming the aircraft at the 
same three airspeeds as for the longitudinal data, but only trim power for level 
flight was used. Four different input forms designated A, B, C, arid D, were 
used to excite the lateral motions, and typical time histories of the inputs 
are illustrated in figure 3, All perturbation data were obtained for one center- 
of-gravity (c,g.) location, 28,8 percent MAC, and test altitudes ranged from 
about 600 m to 1500 m in relatively smooth air. 

The number of test runs made at each condition longitudinally and laterally 
and for each input is presented in table III, 

(2) Steady tests ,- These tests. consisted of trimming the airplane for a 
steady full power climb for a particular airspeed followed by setting the power 
to idle and trimming for an idle power steady descent at the same airspeed 

and altitude range as the climb. This test was done for the same three 
airspeeds as above, with no flap deflection, and two c.g. locations 28,8 percent 
MAC, and 36,5 percent MAC, No analysis for lateral characteristics was 
performed on these data. 
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(3) Quasi-steady tests , » These tests consisted of performing slow 
acceleration-deceleration flight test maneuver (ref. 12), starting from trimmed 
level flight. Power for these tests was left at the trim setting. The initial 
trimmed indicated airspeeds were again 31.6 m/sec (61 knots), 40.7 m/sec 
(78 knots), and 54.2 m/sec (104 knots). These tests were performed at the 
same c.g. locations as the steady tests. No flaps were used nor was an analysis 
of lateral characteristics made for these data. 


DATA REDUCTION AND ESTIMATION METHODS 


The measured flight data included dynamic pressures, angles of attack and 
sideslip, linear accelerations, rotational rates, and control surface movements. 
It was necessary to apply corrections to the measured flight data before it 
could be used for estimating stability and control parameters. Since the angle 
of attack and sideslip vanes were mounted on a boom located near the left wing 
tip and extended about 3/4 chord ahead of the wing leading edge, corrections for 
upwash and angular rates were applied to the measured angle-of-attack. Angle- 
of-sideslip was corrected for angular rates. Details of the angle-of-attack and 
sideslip corrections are given in reference 7. Airspeed was corrected for 
position error by applying a correction to the static pressure determined from 
an airspeed calibration test. The airspeed was also corrected for altitude to 
obtain true airspeed, and since the pitot head was located on the boom, angular 
rates were taken into account to convert airspeed to the aircraft c.g. (ref. 7). 
The accelerometer readings were also corrected to the c.g. of the airplane. 


Three methods were used to estimate stability and control parameters using 
the corrected flight test data. These were the maximum likelihood technique 
described in reference 5, a regression parameter estimation technique described 
in references 3 and 11, and an analytical technique described in reference 13. 


The maximum likelihood technique utilizes the log- likelihood function 


J(0) = - 1/2 E Tl^ r"^ \ - I 


where 
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with z, the measurement vector and y^ the output vector which comes from 
X = f(xt u, e, t) and y = g(x, U, 0, t ) . In the above equation T|^ is 
assumed to have a Guassian distribution and the representation X = f(x, U, 0, t) 
is assumed to accurately represent the physical system. The unknowns to be 
estimated are the elements of 0 and R. Minimizing J with respect to R, 
IT 

R = diag ^ 2 T|^ T|^ is obtained. The estimates for the parameters are obtained 
i 

from the equation 
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yielding the parameter estimates 
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The regression technique utilizes the cost function 
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where r indicates the rth state equation. 

The estimates of the unknown parameters are obtained from the equation 
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se 

which results in 


® - [pi =^ 1 ] *rl] 


where the matrix X. includes measured states and output variables (assumed 
noise free). 


The analytical technique was used to determine the longitudinal static and 
control parameters. These were estimated using the measured stick-fixed trim 
curves at two c.g. positions and the measured as a function of angle-of- 

attack curve. Cx, was determined using the slope of the versus a curve. 


a 


was 


C was determined from the slopes of the stick-fixed trim curves. CL, 
m ^ “‘6 e 

ot 

determined using the difference between the stick-fixed trim curves for various 


Cl’s. 


^6e 


was then calculated from 


(see relation in appendix). 


The application of these techniques is outlined 
likelihood and regression techniques were applied to 
the analytical technique of reference 13 was applied 
quasi-steady data. 


in figure 4. The maximum 
the perturbation data, and 
to the steady and 
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One criterion used to evaluate the uncertainty of the parameter estimates 
obtained using the maximum- likelihood technique is the Craraer-Rao bound 
discussed in reference 14. The Cramer-Rao bounds are an estimate of the 
standard deviations of the parameter estimates but they are too small except for 
the ideal case of a perfect mathematical model, infinite data points, and 
unbiased random noise in the data. Since this is not the case in a practical 
situation, the Cramer-Rao bounds can be used in a relative sense to determine 
the estimation accuracy. 

The effects of power settings on the estimated parameters were 
determined by comparing the values of the estimated parameters determined using 
full power data with those determined using idle power data. 


RESULTS AND DISCUSSION 


Longitudinal 

The results of applying the various estimation techniques to the longi- 
tudinal flight data are shown in figure 5, The longitudinal parameters 
estimated for each longitudinal run are also given in table IV. All of the 
parameters except C^^ and C^^ had Cramer-Rao bounds less than 2 percent of 

the estimated value (table V). These results indicate that the uncertainty in 
the estimated values was "small. Where repeat runs were available, the longi- 
tudinal parameter values estimated by both methods agreed to within 10 percent 
of each other in the majority of the cases. Also, the fit to the flight data 
using the maximum likelihood method was considered good since the mean-squared- 
fit error (area between measured and computed time histories) for each of the 
states was less than 1 percent of the full-scale range of the instrument used to 
measure that state. A typical comparison of measured and predicted flight data 
time histories is shown in figure 6. 

The trends of the estimated parameters with C^ were consistent with those 
approximated from reference 10 where comparisons could be made. The values 
determined by both the maximum likelihood and the linear-regression methods 
generally showed similar trends. The left half of figure 5 shows the parameters 
estimated from data vrtiich were obtained by perturbing the aircraft from trim 
powered level flight while the right half shows parameters estimated from data 
taken when the aircraft was perturbed from idle or full powered fli^t. The 
curves fared through the points determined using the maximum likelihood estima- 
tion method represent the estimates of the derivatives over the Cj^ range for 
which flight data were available. Derivatives determined by other methods were 
shown for comparison, but the values determined using maximum likelihood were 
considered the most reliable. 
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(X 


was positive and increased linearly with This trend seems , 

reasonable since C« was approximately proportional to Ct (ref. 15). There 

a 

was a small power effect at the smallest and middle values tested. 


C_ was approximately constant with Ct, as was expected, since Cx 

varies linearly with a in the a range covered by the flight tests (fig. 5). 
The power effect observed, idle power giving the least negative value, was 
greatest at the largest Cl values. The same trends can be seen in the results 
for both the maximum likelihood and the regression extracted methods. 


Values of for full power, trim power for level flight, and idle power 

were calculated by the techniques of reference 13 using the quasi-steady and 

steady measurements and the values of C„ are shown on figure 5. The values 

cy 

of C^ determined from the steady measurements taken at full power and trim 

power were about the same, and these values of C^^ were more negative than 

those extracted from the ‘perturbation tests. The values calculated by 

the methods of reference 13 from the steady measurements taken at idle power, 
were approximately the same as the values estimated from the perturbation test 
data taken at idle power. 


The power effect on the estimated values of C^^ from steady measurements 

was about 20 percent for this parameter. The power effect on the values 

estimated using maximum likelihood was about 12 percent. As Cl was decreased, 

the difference in the estimated values at C« for full and idle power 

(y 

decreased for both the perturbation and steady measurements. However, the 

values of C« determined from the steady measurements indicated a larger power 

(y 

effect for all C^ values tested. 


Reference 10 describes wind tunnel tests on a high-winged, single-engine 
airplane. Although the configuration tested in the wind tunnel was not exactly 
the same as the configuration of the airplane discussed in this report, the two 
aircraft are somewhat similar. The data in reference 10 indicates that -C^ 




-C^ ) becomes greater as power increases. With the throttle in the 


O' 


<% 

idle position in flight, thrust was -approximately zero, so this case was 

considered to be comparable with the = 0 case of reference 10. The thrust 

coefficient for the full power flight condition was approximately the same as 

the T' = 0.26 case of reference 10. The -C^ of reference 10 for a 
c 


O' 


T* = 0.26 is about 14 percent greater than the -C 2 

C . (y 


value for 


T' = 0. 
c 
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A similar change was noted for the -C^ estimated from flight data using the 

maximum likelihood technique. No reqson has been found for the differences in 
power effects for the values of determined from the perturbation and , 

steady measurements. However, the power effects seen in the results from the 
perturbation tests are of the same magnitude as those seen in the wind tunnel 
test results, so the effect of power on for the maximum likelihood 

estimated results were considered reasonable. 



The magnitudes estimated for and its trend with differed greatly 

for the maximum likelihood and the regression estimation methods. Since the 
parameter had a Cramer-Rao bound at least three times greater than the bound 
for any of the other estimated parameters was not considered well 

determined. Therefore, the estimated value, regardless of its magnitude, did 
not significantly affect the calculated motions of the aircraft. This can be 
seen by examining the period and time to damp to half amplitude of the short 
period mode (table VI). 

While the actual values of all the longitudinal parameters estimated by 
the maximum likelihood ^nd linear regression methods are different, both sets 
give reasonable fits to the flight data. For representative sets of flight data 
the periods and times to damp to half amplitude for the short period mode are 
shown as table VI, The quantities in the table were calculated using the sets 
of parameters determined by examining a particular set of data with both the 
maximum likelihood and regression methods. The difference between the periods 
and times to half were around 10 percent in most cases so that even though the 
individual derivatives in the mathematical model describing the aircraft were 
different the resulting motions were similar. 


The values estimated for C-, using both the maximum likelihood and the 

OS 

regression techniques were similar even though different mathematical models 

were used during the estimation procedure. The maximum likelihood mathematical 

model used a constraint equation which calculated values of Cg from the 

0 s 

estimated value of C_ (see refs. 6 and 7), There were no constraints in the 

*“6e 

regression mathematical model. For the full and idle power cases the magnitudes 
of the estimated parameters were about the same and there were no obvious 
differences in the trends with Cj^ between the two estimation methods. The 
trend with Cj^ seemed somewhat different for the two estimation methods when 
the trim power case was examined, but the magnitudes of the parameters were 
similar (see fig. 5), 
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The estimated values of Cj^ became more negative with increasing Cj^ 

(fig. 5). The trends for Cu^ were the same for both extraction methods, but 

the values determined by the maximum likelihood method were more negative. The 
results showed a definite power effect on Cjj^ for the higher values, 

increasing the power tended to make the values of less negative. 

The estimated values were compared with values calculated from the steady 
and quasi-steady measurements using the methods of reference 13, As can be seen 
from figure 5, the values of Cjj^ calculated from quasi-steady measurements 

showed similar trends and magnitudes as the results of the maximum likelihood 
estimation. 


^mq 

The values of for the full power case tended to be more 

negative as Cj^ increased (fig. 5). Otherwise there was very little variation 
with Cj^. Both parameter estimation methods showed the similar trends and 
magnitudes. The values obtained using the regression method tended to be 
more negative. 

A definite power effect was observed, especially at the largest Cj^, 

This was expected since the dynamic pressure ratio at the tail was greater at 
full power for the largest Cj^ than at full power for the lower Cj^'s, 


The values of ^mg^ determined from perturbation data became more 

negative as increased for both the trim and full power cases and for 

both extraction methods (fig. 5), This increased elevator effectiveness at 
higher Cj^ values is due to a higher d3mamic pressure ratio at the tail for 
trimmed and full power than for the lower Cj^ values. The results from the 
idle power tests showed that remained approximately constant with 

increasing (fig. 5), which would be expected since the propeller slipstream 
is minimal at idle power. The effect of power is similar to that noted for 
^n>q effect of power being greater at the larger Cj^, 

The values determined using the. quasi-steady measurements and maximum 
likelihood estimation showed similar trends with Cj^, For the trim power case 
the Cjjj determined from the quasi-steady measurements showed a much larger 

variation with increasing Cj^ than the Cj^ determined from the maximum 
likelihood estimates, 
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Lateral 


The lateral derivatives estimated by the maximum likelihood method are 
shown on figure 7. A linear fit to each set of these derivatives is also 
shown on figure 7. Values from this linear fit to the derivatives determined 
by maximum likelihood are the preferred values to be used to mathematically 
represent the subject aircraft within the flight regimes covered by the tests. 
Other values are presented for comparison with the values determined using 
maximum likelihood. The lateral derivatives estimated by the linear regression 
method are shown on figure 8, The individual parameters for each run as 
estimated by each method are given as tables VII, and VIII, The run numbers are 
shown on the tables to enable the reader to compare the results of applying 
each method to the same data. For comparison the linear fit to each of the 
maximum likelihood derivatives was also shown on the plot of the corresponding 
derivatives as determined using the linear regression method (see fig, 8). The 
derivative values determined by both methods in general showed similar trends 
and in most cases had similar magnitudes. As another comparison of the maximum 
likelihood and equation error parameter estimates, the characteristics of the 
dutch roll, roll, and spiral modes were computed based on the two sets of 
estimates. These computed characteristics are shown in table IX. The charac- 
teristics estimated are similar for parameters obtained from both techniques. 
The most obvious discrepancy is in the description of the spiral mode, but 
these differences are not considered important since the data runs were not 
long enough to accurately describe the spiral mode. 

Where possible the trends of the parameter values determined using 
maximum likelihood were also compared with trends obtained from references 10 
and 16. Also, the estimated parameter values will be compared with values 
taken from reference 17. These are given in table X. The derivative values 
shown are for a high-winged, single-engine, general aviation aircraft, but only 
the values from reference 17 are for the specific configuration of this report. 


The Cramer-Rao bounds of the lateral derivatives estimated were examined 
and Cyg, C^p, C^.p» ^np» ^nr» *^^6r found to have bounds that 

were less than 2 percent of the extracted values (table XI). An examination of 
figures 7 and 8 reveals that the derivatives listed above had the least scatter, 
as would be expected since the Cramer-Rao bounds were small. In most cases, the 
derivatives Cy^, Cy^^, and had Cramer-Rao bounds of less than 

5 percent of the estimated values, the derivatives C, , and C_ had a 

Cramer-Rao bounds of less than 10 percent of the estimated values, and C„ 

had a Cramer-Rao bound which varied considerably from run to run. The deriva- 
tives with Cramer-Rao bounds of 5 percent or less of the estimated value were 
considered well determined. Also, the fit to the flight data was considered 
good since the mean squared fit error for each state was less than 1 percent of 
the full scale range of the instrument used to measure that state, A typical 
fit to the lateral motions is shown in figure 9, 
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No apparent effect of control input was seen in the values obtained using 
the maxiniuin likelihood estimation method. However, when the linear regression 
estimation method was used, several of the parameters showed some effect of 
input form (see fig. 10). With the exception of Cy^ and none of the 

parameters that were considered well determined showed an effect of input. The 
values estimated for C„a from the data generated when input B was used to 

perturb the aircraft were consistently less negative than when the data generated by 
the other inputs were used. Some of the parameters that were not as well 
determined also had apparent input effects for some C^'s, but the scatter in 
the values estimated by the linear regression method made a definite conclusion 
difficult. 

The reasons for input effects occurring when the linear regression method 
was used and for their absence when using the maximum likelihood method was 
used is unclear. One clue may be in the fact that the cost functions are 
defined differently for the two methods. In the linear regression method, the 
unknown parameters are estimated for each state equation independently of the 
other state equations. So, when using linear regression, since the individual 
inputs excite each of the states in a different manner, several parameter 
values could be estimated to describe the same state. For the maximum likeli- 
hood method all the states to be fitted are estimates simultaneously and the 
unknown parameter values are determined to give a best overall fit to all the 
states simultaneously. 

The parameters Cy^^, and showed good agreement 

between the estimation methods (figs. 7 and 8). These parameters had trends 
with Cl similar to those seen in references 16 and 17 where comparisons could 
be made. With the exception of and the estimated values agreed 

with the values shown in table XI from other references. The lateral parameters 
which had differences between the values determined by maximum likelihood and 
the values determined by the other methods will now be discussed individually. 


C 


ye 



was negative and became less negative as 


Cl increased. 


The maximum 


likelihood results showed no effect of control input and the scatter was small. 
The results using linear regression were less negative than the maximum likeli- 
hood results and the scatter was greater. The general trend for the Cy^ 


values extracted by the linear regression method was the same as for the 
maximum likelihood method, but the linear regression showed some effect of 
input for input B. 
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was not determined well by either estimation method as was indicated 

by the scatter of the estimated values. When using the linear regression 
program, particularly large run to run variations were noted for the runs with 
input B, 


C 




The values estimated for were positive and increased with increasing 

This trend was similar to the trends shown in reference 17, but the magni- 
tudes were lower than expected. The values determined by the linear regression 
showed the same trend as the maximum likelihood but the values estimated for 
the parameters were larger. For both estimation methods the values determined 
using the data for input A were generally lower at the larger Cj^*s than those 
determined using the data from the other inputs. 


C 




The values estimated for were positive and increased with increasing 

Both the maximum likelihood and linear regression estimation methods 
generally resulted in similar trends and magnitudes for the estimated parameter. 
The exception was the values estimated from the data obtained from input A 
using the regression method. These values tended to show a decrease in Cf 
with increasing 


The general trends seen in the estimated values of opposite 

those indicated in reference 16 for a straight-winged aircraft. Also, the 
values estimated were smaller than the value determined from reference 17 and 
shown in table X, However, since is not a strong parameter and is not 

well determined, the apparent discrepancy did not significantly affect the fit 
to the data or the other parameter values estimated. 


The values which were estimated by both methods were negative and 

in general the trend of the derivative values was more negative as 
increased. While both extraction methods gave similar overall trends and 
scatter for the estimated derivatives some effect of input was seen in the 
regression results for input D which showed a reverse trend with and for 

input A which tended to have values less negative than for the other inputs. 

The scatter in the data hid any possible input effect in the maximum likelihood 
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results. When compared to the values given in table X, the magnitudes of Cn 
estimated seemed reasonable and reference 16 Implies that the trend with 
is reasonable. 



The values determined for 


were generally negative, indicating an 

6 3 


apparent proverse yaw with aileron deflection for the sign convention used in 

this report. Experience with the subject aircraft has demonstrated that this 

aircraft actually has adverse yaw with ailerons. This apparent contradiction 

can be explained in part when it is realized that the inputs required to 

determine also produce considerable rolling motion. This rolling motion 

6a 


induces a yawing moment described by the parameter 


Thus, the same 


aileron inputs which cause a yawing moment through also produce a rolling 


motion which causes a yawing moment through C, 


Ilr 


Therefore, the effects 


described by these two parameters is difficult to separate, as is indicated by 

the high correlation (.95) between this pair of parameters. This implies that 

both parameters are required to properly describe the adverse yawing motion 

observed for the subject aircraft, and not or alone. 

“5 a p 


CONCLUDING REMARKS 


The maximum likelihood and a linear regression parameter extraction 
program were used to determine the longitudinal stability and control 
parameters from data taken using a high-winged, general aviation aircraft. 
The parameter values obtained using the maximum likelihood method were 
considered the most reliable representation of the subject aircraft, while 
the values determined by other methods are presented for comparison. 


The longitudinal parameters estimated using maximum likelihood showed the 
trends expected with variations of power and Cl and were considered to give 
a reasonable mathematical representation of the aircraft. For a majority of 
the parameters the same trends were apparent using the linear regression 
method, A definite effect of power setting was observed in the derivatives 


^Z5e» 


+ 


at the largest 


, and 
Cl* The values 


^mgg> with some power effects observed in 


“a 


estimated for C 


'^Zq» -zge* 


and 


C^q + Cj^ were noticeably different when different parameter estimation 


methods were used. The parameter Cj„^^ showed a lesser variation with 

estimation method. Also, the derivatives Cj. , Cg , + Cjj^ , and 

Cm^^ varied with Cj^. -ar o« n- 4 or 
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The parameters ^6e Compared with values determined 

using steady and quasi-steady test data. The trends with were the same 

for these derivatives regardless of the method of determination, but the 
magnitudes were different. The greatest differences were observed in 

where the power effects predicted by the steady and quasi-steady tests were 
different from those predicted by perturbation tests. 


The lateral parameters estimated by both the maximum likelihood and 
equation error methods generally showed the same trends with Cj^, For 


^^p’ ^"p’ ^r’ ^^dr 


the values estimated using the two 


techniques agreed very well. The trends of the parameters estimated were as 
expected and the magnitudes determined were reasonable. As with the longi- 
tudinal results, the parameters estimated using maximum likelihood were 
considered the parameters to be used in any mathematical representation of the 
subject aircraft. The parameters Cn^ and ^n^^ were correlated indicating 

that a description of the adverse yaw due to aileron input of this aircraft 
required a combination of these parameters. 


The estimated parameter values obtained using the maximum likelihood 
extraction method resulted in fit errors less than the uncertainty in the 
measurements for both the longitudinal and lateral data. The estimated longi- 
tudinal parameters all .had Cramer-Rao lower bounds that were less than 


2 percent of their value, as did the lateral parameters 


Cyp, Ci 




"np. 


and Cr 


The lateral parameters C 




and C 


‘ 5 r’ ■ ‘ 'yr* “yer’ “‘'r* “p 

Cramer-Rao lower bounds that were less than 5 percent of the estimated parameter 

had bounds less than 10 percent of the value. 


had 


value, while C 


^6i 


and 


^n 


6a 


Cv had a bound of approximately 25 percent of the extracted value. 
•^P 


The agreement of the two estimation methods, the reasonable trends and 
values of the extracted derivatives, the good fit to the data, and the low 
values of the Cramer-Rao lower bounds gave confidence that for the flight 
data examined, the mathematical model estimated using the maximum likelihood 
method was a reasonable representation of the subject aircraft. 
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APPENDIX 


EQUATIONS OF MOTION 


The equations used in this program are perturbation equations from 
trimmed level flight and are written relative to the set of body axes shown in 
figure 1. 


The equations used to describe the longitudi;nal motions were 

,2, 


u = -qw + rv - g sin 0 + 
w = -pv + qu + g cos 0 cos 




(Al) 


+ (6g - 6e 


66 

(I, - Iv) I 


..>] 


'P + 2 P ^ + % If 


Z X 
q = pr = 

J-xr - 


XZ , 2 2, 


V^Sci 


+ — (r" - P") + P o + - »t) 

L- 


c^c 

+ "^m. 2V ^m 2V ' "m6e"'e 
O' q 


2Iy I m,o “‘a 


0 = q cos q> - r sin co 


(A2) 


(A3) 

(A4) 


a^ = -(u + qw - rv + g sin 0) 

^ 8 

= -^(w + pv - qu - g cos 0 cos cp) 


^u^ + v^ + 


-1 w 
Q- = tan — 


w 


u 


. W 

O' ^ u 


~ Cj„ (used in maximum likelihood extraction) 

^6e 6e 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 
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The values of the lateral states v, p, r, and cp used in the longitudinal 
equations were the flight-measured quantities. 


r ' 


Since thrust changes are not explicitly modeled in the equations of motion. 


O' 


and Cl are not necessarily pure 


and C. 


Of 


contributions due to changes in thrust. Therefore, 
in this study, are given by 


r • = 

^X 

Of 


Sof 


+ C„ 


cos e 


Of 


V but may contain small 

Of 

4 


Of 


and , as determined 

Of 


a (y 


Since, in this study, thrust was held constant and the angle-of-attack changes 
were no more than 7® peak to peak, the contributions of thrust to and 

were considered minimal, cy 

cy 

The equations used to compute the lateral motions were 


V = -ru -f* pw 4* g cos 9 sin 


1 V^sC 
‘P + 2 p ~ S. 


o + P 


+ + 


p 2V 


•r(6r " ^r,t^ 


+ e 


. ^xz . . M . / ^xz\ . 1 V^SbL 

2V ^i-r 2V ^Z'er^^r " ^r,t^ ^z-da^^a " ®a,t^ 

. ^XZ . Ax “ ^ ^ 


o + ‘'np e 


rb 


hp 2V 2V ~ ’’’ 


+ C + 

^ ^np 2V ^ 

9 = p + (q sincp+r cos cp) tan 9 


6a<»a 




a^ = ~(v + ru •- pw - g cos 0 sin cp) 


(All) 


(A12) 


(A13) 
(A 14) 
(A15) 
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+ V + w 


0 = Sin - 


The values of longitudinal states u, w, q, and 0 used in the lateral equa- 
tions were the flight-measured quantities. The equaWons were used to compute 
the airplane state responses. T^e computed ^esponqes were then compared with 
the recorded responses from the flight tests and tlic differences were used to 
update the parameters (stability and control derivatives) to improve the fit. 

The longitudinal measured and computed responses, pr states, used in the 
algorithm for this study were u, w, q, 0, a^^, and a 2 » The lateral states 
used were v, p, r, cp, and ay. Discussion of the identification algorithm 
is given in reference 7, 
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TABLE I.- GEOMETRIC CHARACTERISTICS 


Mass, kg 837.93 

Inertia: 

I , kg-mg 1395 . 

1 ^, kg-nig 11 j 80 . 

I^, kg-m 2 2563 . 

1 x 2 » ^^3. 

Fuselage length, m 8.2 

Wing: 2 

Area, in l6.2 

Aspect ratio T.4T 

Span, m 11.0 

Mean geometric chord, m 1.49 

Vertical tail: 

Area, 1.04 

Aspect ratio 3.96 

Span, m ..... 2.03 

Rudder area, var .68 

Hori zontalgTail : 

Area, m 3.35 

Aspect ratio 3.44 

Span, m 3.45 

Tail length, m 4.36 

a, 3> V Boom location relative to c.g.: 

X, m 1.47 

y> m - 5.43 

m - .77 
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TABLE II.- INSTRUMENT RANGES 


Instrument 

Range 

Airspeed, m/sec 

0 to 63.0 

Angle of attack, deg 

- 8.0 to 39.0 

Angle of sideslip, deg 

, + 23.0 

Altitude, m 

0 to 

Normal acceleration, g units 

-.5 to 4.0 

Longitudinal acceleration, g units 

+ 1.0 

Lateral acceleration, g units 

+ 1.0 

Elevator position, deg 

+25.0 'to -29.0 

Aileron position, deg 

+16.0 to -20.0 

Rudder position, de’g • 

+19.0 

1 

Throttle position 

Total Throttle travel 

Pitch rate, deg/sec 

+30.0 

Roll rate, deg/ sec 

+30.0 

Yav rate, deg/sec 

+ 30.0 

Pitch attitude, deg 

+30.0 

Roll attitude, deg 

+60.0 
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TABLE III.- FLIGHT TEST CONDITIONS 


LONGITUDINAL 


Trimmed 

^“"^^^^Airspeed 

Power 

Setting 

31.6 m/sec 

ltO.7 m/sec 

5^.2 m/sec 

IDLE 

1 run 

2 runs 

1 run 

LEVEL FLIGHT 

2 runs 

2 runs 

2 runs 

FULL 

2 runs 

2 runs 

same runs 
as Level Flight 


LATERAL 


Trimmed 

^^^Airspeed 

Power 

Setting 

31.6 m/sec 

31.6 m/sec 

31.6 m/sec 

31.6 m/sec 

LEVEL FLIGHT 

Input A 
3 Runs 

Input B 
2 Runs 

Input C 
3 Runs 

Input D 
1 Run 

Trimmed 

^\Airspeed 

Power 

Setting 

UO.T m/sec 

1+0. T m/sec 

1 + 0.7 m/sec 

1 + 0.7 m/sec 

LEVEL FLIGHT 

Input A 
1 Run 

Input B 
3 Runs 

Input C 
2 Runs 

1 

Input D 
1 Run 

Trimmed 
^\^Air speed 
Power 
Setting 

5^.2 m/sec 

1 

5 I +.2 m/sec 

5 I +.2 m/sec 

5 I +.2 m/sec 

LEVEL FLIGHT 

Input A 
2 Runs 

Input B 
3 Runs 

Input C 
3 Runs 

Input D 
1 Run 
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TABLE IV.- LONGITUDINAL PARAMETER VALUE 


MAXIMUM LIKELIHOOD ESTIMATION METHOD 
(CRAMER-RAO BOUNDS IN PARENTHESIS) 




-m ^ / r- 










Power Setting 


IDLE 

TRIM 

TRIM 

FULL 

FULL 

Parameter 

c 


1.68 

i. 4 o 

l. 4 i 

1.44 

1.81 

X 

a 


(.0620) 

(.0293) 

(.0183) 

(.0223) 

(. 0247 ) 

c 


-1.076 

-.951 

-.980 

- 1.00 

-.972 

z 

o 


(.002U) 

(.0017) 

( .0030) 

( .0026) 

(.0026) 

p 


-U .28 

- 4,74 

- 4.92 

- 4.76 

- 4.95 

^z 

a 


(. 101 ) 

( .0622) 

(.0577) 

(. 0841 ) 

(.0670) 

c 


- 9.77 

-9.78 

-10.29 

-13.09 

-12.62 

z 

q 


(.969) 

(. 754 ) 

(. 714 ) 

(1.103) 

(.905) 

c 


-.355 

-.425 

-.466 

-.512 

-.496 



- 1.22 

- 1.13 

-1.06 

-.862 

-.947 

m 

a 


(.0206) 

(. 0110 ) 

(.0095) 

(. 0132 ) 

(. 0105 ) 

m 

C ^ 


- 4 . 0 * 

- 4 . 0 * 

- 4 . 0 * 

- 4 . 0 * 

- 4 . 0 * 


-11.82 

- 14.55 

-16.71 

-19.67 

-18. 4 i 

m 

q 


(.298) 

(.188) 

(. 214 ) 

(.272) 

(.218) 

c 


-l.o 4 

-1.25 

- 1.37 

-1.50 

-1.45 

5e 


(. 0113 ) 

(.0071) 

(.0086; 

(. 0102 ) 

(.0083) 


*Parameter Fixed 


LINEAR REGRESSION ESTIMATION METHOD 
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TABLE IV CONTINUED 


MAXIMUM LIKELIHOOD ESTIMATION METHOD 
(CRAMER-RAO BOUNDS IN PARENTHESIS) 



ATPQp-P'P'n Ln T TTi/oor* wl 







^ 

Power 

Setting 

IDLE 

IDLE 

TRIM 

TRIM 

FULL 

FULL 

Parameter 







c 

1.22 

1.01 

.836 

.867 

.937 

.884 

X 

a 

(.0itll2) 

(.0306) 

(.0169) 

(.0208) 

(.0174) 

(.0189) 

c 

-,6Ui 

-.665 

-.645 

-.612 

-.639 

-.6l4 

Z 

0 

(.0026) 

(.0039) 

(.0025) 

(.0022) 

(.0025) 

(.0029) 

c 

-U.6U 

-4.56 

-5.26 

-5.07 

-4.92 

-4.98 

z 

a 

(.0838) 

(. 0947 ) 

(.0472; 

(.0520) 

(.0636) 

(.0722) 

c 

-10. IT 

-10.76 

-8.95 

-9.09 

-11.28 

-10.51 

z 

1 

(1.091) 

(1.149) 

(.585) 

(.642) 

(.855) 

(.929) 

C 

-.361+ 

-.369 

-.422 

-.413 

-.440 

-.440 


-.877 

-.879 

-.853 

-.886 

-.855 

-.825 

■ m 

a 

(.0114) 

(.0155) 

(.0086) 

(.0078) 

( .0098) 

(.0118) 

C 

m. 

a 

C 

-4.0* 

-4.0* 

-4.0* 

-4.0* 

-4.0* 

-4.0* 

-13.00 

-13.51 

-14.94 

-i4.75 

-15.92 

-15.87 

in 

q. 

(.248) 

(.319) 

(.245) 

(.171) 

(.228) 

(.296) 

c 

-1.07 

-1.08 

-1.24 

-1.21 

-1.29 

-1.29 

in ^ 
5e 

(.0090) 

(.0117) 

(.0103) 

(.0062) 

( .0087) 

(.0121) 


*Parameter Fixed 


LINEAR REGRESSION ESTIMATION METHOD 



ATT 3 QPT?T?T^ )\C\ V /c. 






/ « 1 XU / 0 ' 




Power 

Setting 

IDLE 

TRIM 

TRIM 

FULL 

FULL 

Parameter 

c 

•y 

.98 

.85 

.99 

.94 

.90 

a 

C 

2 

-.586 

-.573 

-.56 

-.55 

-.57 

0 

C 

z 

- 4.72 

-4.77 

- 4.72 

- 4.74 

-.485 

cx 

C 

2 

-16.8 

-l 4 .l 8 

-16.27 

-17.11 

-15.39 

q. 

c 

-.36 

-.37 

-.42 

-.46 

-.33 

^6e 

C 

-.66 

-.64 

-.67 

-.63 

VO 

r 

m 

C ^ + c 

-17.76 

-19.0 

-19.5 

-19.76 

- 20.14 

C 

-i.o4 

-1.16 

-1.17 

! 

- 1.21 

-1.23 

6e 
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TABLE IV CONCLUDED 


*Parameter Fixed 


Power 

Setting 


Parameter 


C 


MAXIMUM LIKELIHOOD ESTIMATION METHOD 
(CRAMER-RAO BOUNDS IN PARENTHESIS) 

AIRSPEED ^h ,2 ra/sec 


Power 




Setting 

IDLE 

FULL = 
TRIM 

FULL = 
TRIM 

Parameter 

.652 


.56 

C 

.52 

X 

a 

(.021*1*) 

(.0151) 

( .0221) 

c 

-.339 

-.37 

-.35 

Z 

0 

(.0021*) 

(.0027) 

(.0025) 

c 

-5.02 

-5.18 

-5.27 

Z 

a 

(.0612) 

(.0672) 

(.0553) 

C 

-8.81 

-8.65 

-8.1*8 

z 

q 

(.705) 

(.968) 

(.755) 

c 

-.335 

-.1*05 

-.398 





C 

-.722 

-.68 

-.676 

m 

a 

(.0130) 

(.0101) 

(.0101) 

c 

-i*.o* 

-i*.o* 

-i*.o* 

m. 

a 


-I5.3I+ 


C 

-12.30 

-15.35 

m 

q 

* (.301) 

(.337) 

(.318) 

c 

-.98 

-1.19 

-1.17 

m. 

6e 

( .0112) 

(.0131+) 

(.0127) 


LINEAR REG RESSION ESTIMATION METHOD 
AIRSPEED =5^.2 m/sec 



q or 


.62 
-.357 
-.U.U9 
-19. 5U 
-.38 
-.50 
17.89 

-.985 


FULD= 

TRIM 


-.339 

-U.6U 

-18.87 

-. 1+5 

-.U82 

-19.86 

1.13 


FULL= 

TRIM 


.51+ 

-.336 

-k.62 

-19.56 

-.1*6 

-. 1*85 

- 19 . 1 + 1 + 

- 1.11 
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TABLE V.- THE PERCENT THE CRAMER-RAO BOUND IS OF THE ESTIMATED 
LONGITUDINAL PARAMETER 


Power 

Setting 

Airspeed 

c 

X 

a 

0 

0 

c 

z 

a 

c 

z 

q 

c 

m 

a 

c 

m 

q 

5e 

TRIM 

31.6 m/sec 

3.690 

0.223 

2.360 

9.918 

1.689 

2.521 

1.087 

TRIM 

31.6 m/sec 

2.093 

0.179 

1.312 

7.710 

0.932 

1.292 

0. 568 

TRIM 

U0.7 m/sec 

1.298 

0.306 

1.173 

6.939 

0.896 

1.281 

0.628 

TRIM 

kO.J m/sec 

l.^h9 

0.260 

1.767 

8.426 

1.531 

1.383 

0.680 

TRIM 

5U.2 m/sec 

1.365 

0.267 

1.35»* 

7.171 

1.109 

1.184 

0.641 

TRIM 

5U.2 m/sec 

3.623 

O.U06 

1.806 

10.728 

1.300 

1.908 

o.84i 

FULL 

31.6 m/sec 

3.030 

0.586 

2.077 

10.678 

1.763 

2.347 

1.083 

FULL 

31.6 m/sec 

2.022 

0.388 

0.897 

6.536 

1.008 

1.640 

0.831 

IDLE 

31.6 m/sec 

2.399 

0.359 

1.026 

7.063 

0.880 

1.159 

0.512 

FULL 

U0.7 m/sec 

1.857 

0.391 

1.293 

7.580 

1.146 

1.432 

0.674 

FULL 

ho. 7 m/sec 

2.138 

O.U72 

1.450 

8.839 

1.430 

1.865 

0.938 

IDLE 

U0.7 m/sec 

3.7^^2 

0.708 

1.219 

8.002 

1.801 

2.447 

i.i43 

IDLE 

U0.7 m/sec 

2.903 

0.730 

1.297 

11.191 

1.485 

2.197 

1.126 

IDLE 

5U.2 m/sec 

3. 9*^6 

0.71*t 

1.o49 

8.903 

1.494 

2.072 

1.085 




o 


TABLE VI.- PERIODS AND TIMES TO DAMP TO 1/2 AMPLITUDE 


Flight 

condition 


31.6 m/sec 31.6 m/sec 40.7 m/sec 40.7 m/sec 5^.2 m/sec 54.2 m/sec 
Trim Trim Trim Trim Full Full 

pover power power power power power 


Extraction method M.L. 


Regression M.L. 


Regression M.L. 


Regression 


Period 


1.80 sec 


1.82 sec 1.58. sec 


1.26 sec 1.33 sec 1.94 sec 


Time to damp 
to 1/2 amplitude 


Flight condition 


►295 sec 


,268 sec 



31.6 m/sec 
Idle power 


►235 sec 


.217 sec .174 sec .160 sec 


31.6 m/sec 
Idle power 


54.2 m/sec 
Idle power 


54.2 m/sec 
Idle power 


Extraction method Regression M.L. 


Regression M.L. 


Regression 


Period 


1.90 sec 1.94 sec 1.91 sec 1.29 sec 


Time to damp 
to 1/2 amplitude 


.372 sec .271 sec 


.197 sec 


.173 sec 























TABLE VII.- LATERAL PARAMETER VALUES DETERMINED USING MAXIMUM LIKELIHOOD ESTIMATION 
(CRAMER-RAO BOUNDS IN PARENTHESIS) 


Case 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 ' 

30 

31 

Input 

A 

A 

A 

A 

A 

A 

B 

B 

B 

B 

B 

B 

"l 

1.09 

.975 

1.102 

.633 

.388 

.363 

1.03 

1.10 

.67 

.576 

.562 

.342 

c 

-.51+6 

-.550 

-.51+2 

-.589 

-.595 

-.563 

-.563 

-.536 

-.580 

-.58 

-.57 

-.59 


(.0036) 

(. 001 + 5 ) 

(.0030) 


(.0049) 

(.0042) 

( . 0056) 

(.0060) 

(. 0034 ) 

( . 0052) 

( • 0086) 

( • 0037 ) 

c 

.058 

.060 

.027 

-.089 

-.o 4 i 

-.' 08 l 

.107 

.l 4 o 

.070 

.053 

.042 

.054 


(.0152) 

(.0208) 

(.0155) 


(.0119) 

(.0105) 

(.0172) 

(. 0148 ) 

(.0108) 

( . 0143) 

( . 0288) 

(. 0107 ) 

c 

.150 

. 11+8 

.160 

.147 

.130 

.096 

.152 

.185 

.166 

.127 

.145 

.133 

^Sr 

(. 00 h 6 ) 

(.0071) 

(. 0039 ) 


(.0076) 

(.0052) 

(. 0046 ) 

(. 0045 ) 

(.0054) 

(. 0047 ) 

(.0062) 

(. 0042 ) 

c , 

-.057 

-.059 

-.058 

-.070 

-.071 

-.073 

-.061 

-.062 

-.074 

-.067 

-.061 

-.078 


(.00055) 

(.00051+) 

(.OOO6I+) 


(. 0013 ) 

(.00090) 

( .00067) 

(00071) 

(. 00047 ) 

( .00067) 

(.00089) 

(.00055) 

c. 

-.U05 

-. 1+00 

-.438 

-.430 

-. 4 o 4 

-.44 

-.44 

-.439 

-.484 

-.440 

-.391 

-.477 

P 

(.OOltl) 

(.0036) 

(.0052) 


(.0058) 

(. 0046 ) 

(. 0041 ) 

(. 0039 ) 

(.0033) 

( .oo 4 d) 

(. 0049 ) 

(. 004 l) 

c, 

.0787 

.100 

.128 

.090 

.047 

.079 

.132 

.137 

.108 

.101 

.125 

.078 

ic 

r 

(. 0021 ) 

(.0023) 

(.0023) 


(. 0045 ) 

(.0029) 

(.0022) 

(.0025) 

(. 0020 ) 

(.0022) 

(.0027) 

(. 0019 ) 

c. 

. 00176 

.0051 

.0082 

.0101 

.0030 

.0055 

.0078 

.012 

.0077 

.0058 

.0110 

.0079 

6r 

(.OOOlll) 

(. 0001 + 5 ) 

(. 00045 ) 


(. 0013 ) 

(.00060) 

(. 00041 ) 

(.0039 ) 

(. 00049 ) 

(. 00041 ) 

(.00044) 

(. 00043 ) 

c . 

-, 07 h 2 

-.081 

-.089 

-.096 

-.095 

-.098 

-.101 

-.100 

-.114 

-.104 

-.087 

-.114 

oa 

(.00067) 

(.00065) 

(.00095) 


(. 0013 ) 

(.00099) 

(. 00084 ) 

(.00077) 

(.00069) 

( .00086) 

(.00085) 

( .00088) 

C 

.038 

.038 

.035 

.039 

.044 

.042 

.033 

.038 

.046 

.042 

.031 

.051 

“3 

(. 0001 + 2 ) 

(. 00039 ) 

(. 00045 ) 


(. 00034 ) 

(. 00035 ) 

(. 00033 ) 

(. 00037 ) 

(. 00020 ) 

( . 00034 ) 

(.00052) 

( . 00022 ) 

c 

-.139 

-.IOI+ 

-.124 

-.104 

-.062 

-.071 

-.108 

-.117 

-.058 

-.078 

-.144 

-.031 

n 

P 

(.0029) 

(.0030) 

(. 0040 ) 


(. 0015 ) 

(.0016) 

(. 0021 ) 

(. 0020 ) 

(.0016) 

(. 0020 ) 

(. 0028 ) 

(.0018) 

"nr 

-.121 

(. 0011 +) 

-.117 

(.0015) 

-.109 

(.0016) 

-.101 

-.098 

(. 0013 ) 

-.095 

(. 0010 ) 

-.144 

(. 0012 ) 

-.124 

(. 0013 ) 

-.100 

(.00078) 

-.113 

(. 0010 ) 

-.132 
( . 0019 ) 

-.105 

(. 00082 ) 


-.063 

-.061 

-.062 

-.056 

-.054 

-.051 

-.066 

-.061 

-.057 

-.059 

-.060 

-.057 

5r 

(.00031+) 

(. 00037 ) 

(.00032) 


(. 00037 ) 

(.00032) 

(. 00035 ) 

(.00036) 

(.000 28 ) 

(.000 30 ) 

(.00043) 

(.000 24) 

C 

-.011 

-.0082 

-.012 

-.012 

-.0061 

-.0068 

-.0015 

-.0084 

-.0011 

-.0029 

-.0134 

.0019 

“6a 

(.00059) 

(.00058) 

(.00078) 


(. 00033 ) 

(. 00037 ) 

(. 00048 ) 

(. 00044 ) 

(.000 38 ) 

(. 00042 ) 

(. 00042 ) 

(. 00041 ) 


TABLE VII.- CONTINUED 


Case 

32 

33 

34 

35 

36 

37 

38 

39 

4 o 

4 i 

42 

43 

44 

Input 

B 

B 

C 

c 

C 

C 

C 

c 

c 

c 

c 

C 

c 

"l 

.320 

.325 

.326 

• 322 

.326 

.385 

.376 

.568 

.662 

.949 

1.00 

.988 

.926 

c 

-.58 

-.58 

-.582 

-.583 

-.581 

-.577 

-.577 

-.602 

-.596 

-.592 

-.59 

-.574 

-.53 


(.0037) 

(.0036) 

(. 0037 ) 

(. 0035 ) 

(.0036) 

(.0033) 

( .0030) 

( .oo 4 o) 

( .0032) 

(. 0054 ) 

( .0032) 

(.0026) 

(. 0044 ) 

c 

.0057 

.o 4 i 

-.122 

-.092 

-.116 

-.025 

-.011 

-.181 

-.122 

-•OTO 

-.079 

-.032 

.210 


(. 0148 ) 

(. 0154 ) 

(. 0102 ) 

(.0098) 

( .0100) 

( . 0110 ) 

(. 0113 ) 

(.0139) 

(.0119) 

(.0195) 

(.0132) 

(. 0104 ) 

( .oi 4 i) 

c 

.lOU 

.121 

.091 

.111 

.095 

.114 

.121 

.134 

.125 

.178 

.155 

.139 

.240 

^6r 

(. 0057 )' 

(.0067) 

(. 0054 ) 

(. 0054 ) 

(. 0054 ) 

(.0051) 

( . 0044 ) 

(.0055) 

( .0050) 

(. 0064 ) 

(. 0042 ) 

(.0039) 

(.0075) 

c. 

-.087 

-.081 

-.074 

-.074 

-.072 

-.081 

-.082 

-.069 

-.0705 

-.053 

-.059 

-.064 

-.064 

3 

( .00088) 

(.00088) 

(.00066) 

(.00072) 

(.00055) 

( .00074) 

(.00091) 

( .00065) 

( .00069) 

( .00066) 

(. 00042 ) 

(. 00041 ) 

(.00059) 


-.557 

-.505 

-.461 

-.457 

-.456 

-.49 

-.485 

-.448 

-.442 

-.377 

-.416 

-.445 

-.45 

36 

p 

(.0072) 

(. 00 t 2 ) 

( .oo 4 i) 

(. 0044 ) 

(. 0034 ) 

(. 0048 ) 

(.0059) 

(. 0045 ) 

( . 0048 ) 

(. 0048 ) 

( . 0048 ) 

( .0036) 

( .0032) 

C, 

.0U7 

. 059 ' 

.076 

.080 

.071 

.062 

.078 

.111 

.093 

.161 

.124 

.113 

.128 

z 

r 

(.0029) 

(. 0037 ) 

(.0021) 

(.0026) 

(.0018) 

(. 0022 ) 

(.0025) 

( . 0020 ) 

( .0026) 

(.0028) 

(.0016) 

(.0016) 

( .0026) 


-.0005 

.0045 

.0035 

.0060 

.0023 

.0036 

.0094 

.0083 

.0087 

.013 

.0094 

.0071 

.0125 

36 . 

6r 

(.00068) 

(.00090) 

(. 00049 ) 

(.00058) 

(. 00042 ) 

(.00055) 

(.00057) 

( . 00049 ) 

(.00062) 

( .00065) 

(.00045) 

(. 00042 ) 

(.00057) 


-.127 

-.117 

-.103 

-.105 

-.102 

-.114 

-.115 

-.102 

-.103 

-.094 

-.099 ’ 

-.103 

-.114 

oa 

(.0015) 

(.0015) 

(. 00084 ) 

(.00092) 

(.00069) 

(. 0010 ) 

(.0013) 

(.00091) 

(. 0010 ) 

(. 0010 ) 

(. 00077)1 

(.00065) 

(.00079) 

c 

.051 

.051 

.044 

.044 

.045 

.049 

.046 

.040 

.o 4 o 4 

.034 

.038 

.029 

.034 


(. 00030 ) 

(.00027) 

(. 00022 ) 

(.00023) 

(.00019) 

(.00023) 

(.00027) 

( . 00030 ) 

( .00029) 

(. 00041 ) 

( . 00022 ) 

( . 00024 ) 

( .00026) 

c 

-.028 

-.030 

-.063 

-.053 

-.061 

-.037 

-.055 . 

-.090 

-.086 

-.095 

-.067 

-.130 

-.073 

n 

P 

(.0026) 

(.0025) 

(. 0014 ) 

(. 0014 ) 

(. 0012 ) 

(.0016) 

(.0019) 

( . 0023 ) 

(.0023) 

(. 0034 ) 

(.0019) 

(. 0020 ) 

(.0019) 

Cti 

-.101 

-. 111 ' 

-.096 

-.093 

-.096 

-.100 

-.101 

-.099 

-.104 

-.125 

-.137 

-.139 

-.094 


(. 0010 ) 

(. 0012 ) 

(.00069) 

(.00075) 

(.00063) 

(. 00064 ) 

(. 00064 ) 

( .00089) 

(.0010) 

(.0016) 

(.00083) 

(.00090) 

( . 0011 ) 

c 

-.055 

-.057 

-.052 

-.051 

-.054 

-.054 

-.053 

-.056 

-.056 

-.062 

-.064 

-.062 

-.056 

or 

(.00028) 

(.00032) 

(.00029) 

(.00029) 

(.00026) 

(. 00023 ) 

(.00021) 

( .00032) 

(.00031) 

(. 00045 ) 

(.00026) 

(.00025), 

( .00027) 

c 

.0014 

.0029 

-.0054 

-.0037 

-.005 

+.0013 

-.0035 

-.0081 

-.0075 

-.00036 

.0069 

-.0064 

.0015 

^6a 

(.00059) 

(.00057) 

(.00030) 

(.00031) 

(.00026) 

(.00036) 

(. 00043 ) 

(.00053) 

(.00051) 

(.00077) 

(. 00045 ) 

(. 00044 ) 
' 1 

( . 00045 ) 


I 




TABLE VII.- CONCLUDED 


Case 

U 5 

56 

Input 

D 

D 

"l 

.559 

.320 

C 

-.56 

(.00U6) 

-.59 

(.0032) 

c 

.055 

(.01U3) 

-, 04 l 

(. 0084 ) 

c ■ 

^ 6 r 

.13U 

(.ootH) 

.131 

(. 004 i) 

3 

-.076 

(.00095) 

-.078 

(.00073) 

P 

-.hj 

(.0058) 

-.47 

(. 0042 ) 

r 

.080 

(.0038) 

.071 

(.0025) 

^ 6 r 

.0114 

(.00062) 

.0097 

(. 00048 ) 

6 a 

-.113 

(. 0012 ) 

-.108 

(.00090) 

C 

“3 

.OU3 

(. 00033 ) 

.047 

(. 00021 ) 

C 

n 

P 

-.072 

(.OO2U) 

-.049 

(. 0013 ) 

Cn, 

-.124 

(. 0015 ) 

-.118 

(.00079) 

S 

5 r 

-.065 

(.00038) 

-.058 

(. 00021 ) 

C 

6 a 

.00036 

(.00052) 

-.0015 

(.00028) 




CO 

4 ^ 


TABLE VIII.- LATERAL PARAMETER VALUES DETERMINED USING 
LINEAR REGRESSION ESTIMATION 


1.102 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

A 

A 

A 

B 

B 

B 

B 

B 

B 

B 

B 

.663 

.388 

.363 

1.03 

1.10 

.67 

.576 

.562 

.342 

.32 

.325 

-.1+76 

-.533 

-.534 

-.388 

-.322 

-.352 

-.373 

-.303 

-.437 

-.432 

-.42 

-.012 

.033 

.0202 

.214 

.349 

.0864 

.377 

.666 

.100 

.034 

.023 

.119 

.125 

.127 

.138 

.137 

.109 

.0908 

.102 

.035 

.011 

.030 

-.0608 

-.0645 

-.066 

-.0476 

-.047 

-.0603 

-.0637 

-.054 

-.075 

-.077 

-.075 

-.378 

-.370 

-.388 

-.365 

-.357 

-.383 

-. 4 o 8 

-.370 

-.44 

-.45 

-.45 

.076 

.0509 

.045 

.133 

.129 

.082 

.115 

.149 

.077 

.067 

.076 

-.0873 

-.088 

-.091 

-.0902 

-.087 

-.091 

-.10 

-.094 

-.108 

-.106 

-.106 

.0049 

.0031 

.00466 

.0093 

.0087 

.0043 

.010 


.0095 

.0074 

.0087 

.036 

.0401 

.0426 

.026 

.035 

.036 

.o 4 o 

.047 

.048 

.048 

.048 

..0708 

-.0485 

-.038 

-.104 

-.0745 

-.0674 

-.080 

-.050 

-.052 

-.051 

-.057 

-.098 

-.092 

-.092 

-.125 

-.102 

-.0985 

-.127 

-.110 

-.117 

-.119 

-.123 

..0062 

-.00309 

-.00068 

-.005 

i 

-.0013 

-.0046 

-.0044 

-.0015 

-.0034 

-.0036 

-.0036 

-.0512 

-.0488 

-.0477 

-.057 

-.053 

-.050 

-.061 

-.058 

-.058 

-.058 

-.059 















































Case 

3 h 

35 

36 

37 

38 

39 

Input 

C 

C 

c 

c 

C 

c 

"l 

.326 

.322 

.326 

.385 

.376 

.568 

c 

-.508 

-.516 

-.509 

-.502 

-.515 

-.500 








c 

-.056 

-.068 

-.073 

-.028 

-.013 

-.065 








p 







c 

.066 

.073 

.069 

.082 

.087 

.124 

^ 6 r 





C, 

-.076 

-.075 

-.076 

-.076 

-.077 

-.071 








Cji 

-.lt 4 ■ 

-.45 

-.45 

-.45 

-.45 

-.453 

p 

.069 





.096 


.072 

.070 

.071 

.073 

r 








-.106 

-.106 

-.106 

-.108 

-.109 

-.108 

^ 6 a 






6 r 

.0078 

.0073 

.0075 

.0071 

.0087 

.0064 

C 

. 01^6 

.046 

.047 

.046 

.045 

.042 

“e 






c 

Y 1 

-.054 

-.049 

-.0503 

-.'059 

-.062 

-.069 

XI 

p 







c 

n 

-. 114 - 

-.116 

-.115 

-.12 

-.118 

-.118 

r 

C 

-.0031 

-. 0021 

-.0025 

-.0033 

-.0045 

-.0034 

“ 6 a 







C 

-.055 

-.056 

-.056 

-.057 

-.056 

-.058 

n, 

6 r 








4 o 

4 i 

42 

43 

44 

45 

46 

c 

c 

c 

c 

D 

D 

D 

.662 

.949 

1.00 

.988 

.926 

.559 

.320 

..49 

-.447 

-.452 

.464 

-.537 

-.56 

-.552 

.061 

.085 

.012 

.081 

.047 

-.037 

-.057 

.118 

.155 

.159 

.144 

.14 

.082 

.060 

■ .071 

-.057 

-.058 

-.064 

-.060 

-.073 

-.083 

-.451 

-.412 

-.421 

-.438 

-.44 

-.475 

-.48 

.092 

.145 

.134 

.117 

.134 

.094 

.074 

00 

0 

H 

-.102 

-.104 

-.105 

-.111 

-.115 

-.•ll 4 

.0068 

.0102 

.0085 

.0047 

.012 

.0075 

.0103 

C\J 

0 

.034 

.035 

.029 

.0302 

.042 

.047 

•.0709 

-.094 

-.086 

-.126 

-.101 

-.078 

-.056 

•.119 

-.138 

-.14 

-.143 

-.111 

-.122 

-.127 

•. 004 l 

-.0021 

-.'000046 

-.0068 

-.0045 

-.0028 

-.0029 

•.058 

-.066 

-.067 

-.063 

-.060 

-.065 

-.061 
















































TABLE IX.- CHARACTERISTICS CALCULATED FROM ESTIMATED LATERAL PARAMETERS 


MODE 

^L 

ESTIMATION METHOD 

PERIOD (SEC) 

TIME TO HALF OR DOUBLE 
AMPLITUDE (SEC) 

Dutch Roll 

1.0 

Maximum Likelihood 

4.263 

1.495 


1.0 

Linear Regression 

4.370 

1.600 


.6 

Maximum Likelihood, 

3.402 

1.497 


- 6 

Linear Regression , 

3 . 4 i 6 

1.522 


.35 

Maximum Likelihood 

2.598 

1.337 


.35 

Linear Regression 

2.611 

1.382 

Roll 

1.0 

Maximum Likelihood 

— 

.162 


1.0 

Linear Regression 


.174 

• 

.6 

Maximum Likelihood 


.122 


.6 

Linear Regression 


.126 


.35 

Maximtun Likelihood 


.092 


.35 

Linear Regression 


.092 

Spiral 

1.0 

Maximum Likelihood 


43.95* 


1.0 

Linear Regression 


57 . 98* 


.6 

Maximum Likelihood 


81.27 


.6 

Linear Regression 


254.85 


.35 

Maximum Likelihood 


4 i .70 


.35 

Linear Regression 


46.00 


*Time to douBle amplitude 




TABLE X.- COMPARISON OF DERIVATIVE VALUES DETERMINED BY THE MAXIMUM LIKELIHOOD 
METHOD WITH THOSE FROM SELECTED REFERENCES 



Reference 

10 

Reference 

17 

Reference 

16 

Values from least squares' linear 
fit to the M.L. estimated para- 
meters for three C values 

L 


Cl = 1.1 

Cl = .25 

Cl= .9 

^L= -3 


Cl = 1.0 

,c 

- . 5 h 

- .22 

- -.303 

- .59 

- .57 

- .55 

C 


- .01+7 

- .213 

- .0I5 

- .005 

.01 

c 

.lU 3 

.15 


.11 

.ll 

.17 

^8 

- .0785 

- .082 

- .122 

- .078 

- .070 

- .058 

P 


- .17 

- .I9I 

- .17 

- .15 

- .ll 

r 


.07 


.07 

.095 

.13 

5 r 


.013 


.0058 

. 007 

.0095 


- . 102 

- .198 

I 

- .11 

- .102 

- .094 

C 

^8 

.057 

.035 

. 0701 

.018 

.0I2 

.033 

C 

n 

p 


- .077 

- .096 

- .05 

- .078 

- .115 

C 

n- 

r 


- .062 

- .115 

- .10 

- .11 

- .125 

C 

n. 

6 r 

- .057 

- .016 


- .053 

- 0.58 

- .063 

C 

n* 

5 a 


. oil 


- .003 

- .001 

- .005 






CASE NUMBER 


PARAMETER 

20 

C 

0.659 

c 

26.21 

c 

3.220 

C 

3.067 


0.965 

3 

•1.012 

P 

2.668 

r 

C, 

23.30 

^<Sr 

0.903 

6a 

C 

1.105 

“3 

C 

2.086 

n 

P 

C 

.1.157 

n 

r 

C 

0.5^0 

“6r 

^“6a 

5.36U 


21 22 


0.818 

0.554 

34.67 

57.41 

4.692 

2.507 

4.798 

2.438 

0.915 

1.103 

0.900 

1.187 

2.300 

1.797 

8.824 ' 

5.488 

0.802 

1.067 

1.026 

1.286 

2.885 

3.226 

1.282 

1.468 

0.607 

0.516 

7.073 

6.500 


0.824 

0.746 

29.02 

12.96 

5.800 

3.937 

5.846 

5.417 

1.831 

1.233 

1.436 

1.045 

9.574 

3.671 

43.33 

10.91 

1.368 

1.010 

0.773 i 

0.833 

2.4 i 9 ' 

2.254 

1.327 , 

1.053 

0.685 ’ 

0.627 

5.4 iO 

5.44 i 


26 


27 


28 


29 


0.995 

16.07 

U.305 

3.026 

1.098 

0.932 

1.667 

5.256 

0.832 

1.000 

1.9^^ 

0.833 

0.530 

32.000 


1.12 

0.586 

10.57 

15.43 

3.963 

3.00 

2.432 

3.25 

1.145 

0.635 

0.888 ^ 

0.682 

1.825 

1.852 

3.25 

6.364 

0.770 

0.605 

0.974 

0.435 

1.709 

2.759 

1.048 

0.780 

0.590 

0.491 

5.238 

34.55 


0.897 

26.98 

3.890 

3.701 

1.00 

0.909 

2.178 

7.069 

0.827 

0.810 

2.56U 

0.885 

0.508 

IU.U8 






PARAMETER 


30 


31 


32 


33 


c 

1.590 

0.508 

0.638 

0.621 

0 

68.57 

19.81 

259.65 

37.56 

c 

5.599 

2.808 

4.838 

4.452 


U.2T6 

3.158 

5.481 

5.537 


1.1*59 

0.705 

1.011 

1.086 

"il 

P 

1.253 

0.860 

1.293 

1.426 

r 

2.160 

2.U36 

6.170 

6.271 

Cji 

^6r 

U.OOO 

5. 1*1*3 

136 

20.0 

*6 a 

0.977 

0.772 

1.181 

1.282 

C 

“b 

1.677 

0.1*31 

0.588 

0.529 

C 

n 

P 

1.91*1* 

5.806 

9.286 

8.333 

C 

n 

r 

1.1*39 

0.781 

0.990 

1.081 

Sr 

‘0.717 

0.421 

0.509 

0.561 

6a 

3.13I* 

21.58 

42.14 

19.66 


CASE NUMBER 

34 35 

36 

37 

38 

0.636 

0.600 

0.620 

0.572 

0.520 

8.36 

10.65 

8.621 

44.00 

102.7 

3.69 

3.50 

3.589 

3.573 

2.911 

5.93 

4.775 

5.579 

4.474 

3.636 

0.892 

0.973 

0.764 

0.914 

1.110 

0.889 

0.963 

0.746 

0.980 

1.196 

2.763 

3.250 

2.535 

3.548' 

3.205 

14.00 

9.667 

18.261 

15.28 

6.064 

0.816 

0.876 

0.676 

0.877 

1.130 

0.500 

0.523 

0.422 

o'.k69 

0.587 

2.222 

2.642 

1.967 

4.324 

3.455 

0.719 

O'. 806 

0.656 

0.640 

0.634 

0.538 

0.549 

0.481 

0.426 

0.396 

5.556 

8.378 

5.200 

27.69 

12.29 






o 


TABLE XI - CONCLUDED 


PARAMETER 

39 

C 

0 . 66 k 

C 

7.680 

c 

3.58U 

C 

U.ioU 

y6r 

C„ 

0.9^t2 

^3 

c„ 

1.004 

P 

C„ 

■ 1.802 

r 

5.904 


0.892 

6a 

C 

0.T50 

“3 

C 

2.556 

n 

P 

C 

0.899 

n 

• 0.5T1 

^“6a 

6.543 


40 

4 l 

0.536 

0.912 

9 -T 54 

27.86 

3.972 

4.337 

4.000 

3.596 

0.979 

1.245 

1.086 

1.273 

2.796 

1.739 

7.126 

5.000 

0.971 

1.064 

0.718 

1.206 

2.674 

3.579 

0.962 

1.280 

0.554 

0.726 

6.800 

213.89 


CASE NUMBER 
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Figiire k.- Flight data processing scheme. 
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Figure 10. - Effect of control input form on estimated parameter values for the lateral inputs shown in 

figure 3. 
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Figure 10. - Continued, 
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Figure 10.- Concluded. 


©• o» 



2. Government Accession No. 


3. Recipient's Catalog No. 


1. Report No. 

NASA TM-80163 


4. Title and Subtitle 

COMPARISON OF STABILITY AND CONTROL PARAMETERS FOR A 

5. Report Date 

September 1979 

LIGHT, SINGLE -ENGINE, HIGH-WINGED AIRCRAFT USING 
DIFFERENT FLIGHT TEST AND PARAMETER ESTIMATION TECHNIQUEi 

6. Performing Organization Code 

7. Author(s) 

William T, Suit and Robert L, Cannaday 

8. Performing Organization Report No. 


10. Work Unit No. . 

9. Performing Organization Name and Address 

505-06-63-02 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

' 

Technical Memorandum 

14. Army Project No. 

15. Supplementary Notes 

16. Abstract 


Longitudinal and lateral .stability and control parameters were estimated from 
flight data for a high-wing, general aviation, airplane using flight data obtained 
at various flight conditions within the normal range of the aircraft. These 
parameters were estimated using an output error technique (maximum likelihood) and 
equation error technique (linear regression). Longitudinal static parameters were 
also estimated from climbing, descending, and quasi-steady-state flight data. For 
the lateral excitations, four input forms were used involving some combination of 
rudder and ailerons. The resulting longitudinal and lateral parameter estimates 
were used to compute the periods and time-to-damp to one-half amplitude of the 
various aircraft modes of motion to determine the sensitivity of these motions to 
variations in the parameter estimates. 


17. Key Words (Suggested by Author(s)) 18, Distribution Statement 

Parameter estimation 

Stability and control parameters Unclassified - Unlimited 

Maximum likelihood 

Subject Category 08 

19. Security Oassif. (of this report) 20. Security Classif. (of this page) 21. No. of Pages 22. Pfice 

Unclassified Unclassified 52 $5.25 

For sale by the National Technical Information Service, Springfield, Virginia 221 61 















I 





